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ABSTRACT 

The disruption of the M33 galaxy is evident from its extended gaseous structure. We present new 
data from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) Survey that show the full extent 
and detailed spatial and kinematic structure of M33's neutral hydrogen. Over 18% of the HI mass of 
M33 (Mnitot = 1.4 X 10^ Mq) is found beyond the star forming disk as traced in the far-ultraviolet 
(FUV). The most distinct features are extended warps, an arc from the northern warp to the disk, 
diffuse gas surrounding the galaxy, and a southern cloud with a filament back to the galaxy. The 
features extend out to 22 kpc from the galaxy center (18 kpc from the edge of the FUV disk) and the 
gas is directly connected to M33's disk. Only five discrete clouds (i.e., gas not directly connected to 
M33 in position-velocity space) are catalogued in the vicinity of M33, and these clouds show similar 
properties to Galactic and M31 halo clouds. M33's gaseous features most likely originate from the 
tidal disruption of M33 by M31 1-3 Gyr ago as shown from an orbit analysis which results in a 
tidal radius < 15 kpc in the majority of M33's possible orbits. M33 is now beyond the disruptive 
gravitational influence of M31 and the gas appears to be returning to M33's disk and redistributing 
its star formation fuel. M33's high mean velocity dispersion in the disk (^18.5 km s^'^) may also 
be consistent with the previous interaction and high rate of star formation. M33 will either exhaust 
its star formation fuel in the next few Gyrs or eventually become star formation fuel for M31. The 
latter represents the accretion of a large gaseous satellite by a spiral galaxy, similar to the Magellanic 
Clouds' relationship to the Galaxy. 

Subject headings: galaxies: individual (M33) — galaxies: halos — Local Group — galaxies: kinematics 
and dynamics — galaxies: ISM — galaxies: evolution 



1. INTRODUCTION 

Large spiral galaxies at low redshift are thought to 
gradually evolve through the accretion of material from 
their halos. This halo material originates from the ac- 
cretion of satellite galaxies and the cooling of the hot 
gaseous halo, with the hot halo most likely originating 
from a combination of the initial baryon collapse and 
galaxy feedback. Continual gas accretion provides star 
formation fuel and allows for the range of stellar popu- 
lations and metallicities we find in galaxies. Simulations 
that examine the methods galaxies obtain gas indicate 
the dominant mechanism depends on a galaxy' s total 
mass and redshift ( Keres et al . 2009, 2005; B rooks et al.l 
[20f)9h . At z = massive galaxies are expected to accrete 
gas through satellites and the gradual accretion of hot 
halo gas, while less massive galaxies (< few x lO'^^ Mq) 
obtain star formation fuel primarily via direct accretion 
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of cold gas along cosmic filaments. 

Observational evidence for the various accretion modes 
remains limited. The Galaxy and M31 have been 
shown to have numerous halo clouds (e.g., MuUer et 
al. 1963; Putman et al. 2002; Thilker et al. 2004), 
with their origin definitely partly from satellite accre- 
tion (|Mathewson et al.lll97i iPutman et al.ll2003D . and 
possibly also from condensing halo clouds in the ho t 
halo medium (!Peek et al.l 12008: Mailer & Bullock 2004). 
There has been additional evidence for halo clouds that 
will fuel spiral g alaxies through deep HI observations 
of ne arby spirals (jOosterloo et al.|[2007l : iFraternali et al.l 
120011 ). In all of these cases, the amount of fuel cur- 
rently found in the halo is not large considering the 
expected star formation rates and chemical evolution 
models. For instance, the Milk y Way is expected to 
have an average i nfall o f 1 Mo/yr (|ChiaPDini et al.l[200ll : 
iFenner fc GibsonI l2003f) . but this reservoir is not read- 
ily apparent from either the accretion of gas from small 
satellites (Grccvich & Putman 2009) or the condensing 
halo clouds (jPeek et al.l I2008f ). In the case of galaxies 
smaller than the Milky Way, observational evidence for 
halo clouds remains limited, and it is possible they will 
largely not be found, as the bulk of the accretion may 
have been completed for these systems at early times. 

M33 is often considered a dwarf spiral galaxy of the Lo- 
cal Group, with a total mass of ~ 5 x 10^" M0(Corbelli 
2003), which is ~20 times lower than the Milky Way 
and M31, and only ~2 times higher than the com- 
bined mass of the Magellanic Clouds. M33 is ~ 730 
kpc from the Milky W ay and ~ 200 kpc from M31 
(jBrunthaler et al.ll2005l ). Its stellar component indicates 
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it is a fairly quiescent galaxy with no evidence for a 
recent interaction. Deep optical surveys have found a 
stellar halo component, but no clear evidence for the 
stream-like sub-structures found in the halo of the Milky 
Way and M31 (Ferguson et al. 2007; McConnachie et al. 
2006). In neutral hydrogen there is evidence for a strong 
warp that has been examined by iCorbelli Sz Schneideii 
(|1997[ ) (hereafter CS97) and Rogstad, Wright & Lock- 
hart (1976). Both authors concluded that the tidal 
force from M31 is unlikely to be sufficient to generate 
the warp. On the other hand, recent work estimating 
the motion of M31 from indirect methods indicates that 
M33 and M31 are almost certainly gravitationally bound 
and M33 should show signatures of tidal perterbatio n 
(jvan der Marel fc Guhathakurtall2008l : iLoeb et al.|[2005D . 

Sensitive large area neutral hydrogen surveys of M33's 
extended halo can investigate the methods that M33 is 
losing and/or obtaining star formation fuel, and the fate 
of this galaxy in the context of galaxy mergers at z = 0. 
Recent observations taken as part of the GALFA-HI Sur- 
vey provide these data at unprecedented resolution (3.4', 
0.18 km s~^) and sensitivity (0.1 K (3ct) over 5 km s^^).^ 
See also Grossi et al. (2008) for data taken with a differ- 
ent spectrometer and reduced with different techniques 
as part of the ALFALFA survey. In this paper we present 
the GALFA-HI Survey observations of the HI environ- 
ment of M33, complete simulations of the orbit of M33 
relative to M31, and discuss the results in the context of 
fueling both M33 and M31. The gaseous environment of 
M33 is an important tracer of the past and future evolu- 
tion of the galaxies of the Local Group. 

2. OBSERVATIONS AND DATA REDUCTION 

The HI data presented here were obtained with the 
Galactic Arecibo L-Band Feed Array (GALFA) spec- 
trometer on the Arecibo Radio Telescope. The Arecibo 
L-Band Feed Array (ALFA) is a 7-beam array of receivers 
similar in design to the 13-beam array used for the HI 
Parkes All-Sky Survey (HIPASS; Barnes et al. 2001). 
GALFA data provide a channel width of 0.18 km s^^ and 
cover a velocity range of —753 to -1-753 km s^^ in the 
Local Standard of Rest (LSR) frame. The spatial resolu- 
tion is approximately 3.4'. At the distance to M33 (730 
kpc), 3.4' corresponds to 720 parsecs. The GALFA-HI 
Survey will eventually map the Arecibo sky from Deck 
= — 1 — 38° and the full RA range. The data presented 
here are from the TOGS (Turn On GALFA Survey) por- 
tion of the GALFA-HI Survey. TOGS obtains data by 
commensally observing with the extragalactic HI surveys 
(ALFALFA (e.g., Giovanelh et al. 2005) and AGES (e.g., 
Auld et al. 2006)), with the majority of the data shown 
here taken commensally with ALFALFA. The observing 
method used is drift scans with our least-squares fre- 
quency switching calibration procedure run at the start 
of each observing run. There are also some basketweave 
scans (driving the telescope only in zenith) taken as part 
of the A2050 GALFA proposal (Pi's: Peek & Heiles) in- 
tegrated into the final data cube that improve any drift 
scan artifacts and increase the sensitivity. 

The reduction process and data are discussed in Peek & 
Heiles (2008) and include corrections for the IF bandpass 

^ Th e data are publicly available through a link at 
|http: / /sites. google.com/site/galfahi/ 



and gain variation in the receivers. The dataset shown 
here includes a sidelobe correction. ALFA is known to 
have significant first sidelobes, which have been mapped 
and are described in Heiles (2004). These sidelobes are 
displaced about 5' from the main beams and have a typ- 
ical sidelobe efficiency of 13% for the central beam and 
23% for each of the six exterior beams. To correct for 
this effect, we use an algorithm similar to the one used in 
Hartmann et al. (1996) for the Leiden-Dwingeloo Survey 
(LDS). We cannot use that algorithm in its exact form, 
as we sample the sky with 7 different beams rather than 
only one beam as in the LDS. In our modified method 
we determine the contamination to each of our observed 
spectra from the first sidelobe. We do this by deter- 
mining the position each sidelobe has on the sky during 
each integration and calculating the spectrum the side- 
lobe would have added to the main beam data. We then 
subtract the 'sidelobe contamination spectra' from orig- 
inal spectra, scaled by the beam efficiency, and use the 
'decontaminated' spectra to construct a new map. 

The data we started with for examining M33 and 
its surroundings were gridded into a cube smoothed to 
0.74 km s^^ velocity resolution. To bring out diffuse 
halo features we also created another cube smoothed to 
5.15 km resolution. This latter cube was used for 
the majority of the analysis presented here and has a 
3(7 sensitivity of 0.1 K per channel, or 10^^ cm~^ in 
each ~ 5 km s""'^ channel. The typical linewidth of 
halo clouds of the Milky Way is approximately 25 km 
s-i (|de Heii et all [2002 ). indicating we are sensitive to 
these types of clouds with column densities above 5x10^^ 
cm~^ {3a-). The typical peak column density of Galac- 
tic HVCs is - 10^9 cm-2 (jPutman et al.|[200l ). The 
3(7 mass sensitivity to objects with linewidths of 25 km 
s~^ at 730 kpc is 2.4 x 10^ Mq. This mass is generally 
at least an order of magnitude lower than the Galactic 
HVC complexes with distance constraints (i.e., Thom et 
al. 2006, 2008; Wakker 200 1) and the lowest HI mas s 
Local Group dwarf galaxies ijGrcevich Sz Putmanll2009l ). 
Before making the figures that include the entire velocity 
range of M33, the emission from M33 was isolated in the 
south where it begins to overlap with Galactic emission 
(see Figure [3]). Diffuse M33 features to the south of M33 
that are indistiguishable from Galactic emission may be 
lost; however the extremely high velocity resolution of 
the GALFA-HI data make it the best data available for 
distinguishing M33 features from Galactic emission. 

3. RESULTS 

3.1. GALFA-HI Observations 

The integrated intensity HI map of M33 is shown in 
Figure [1] with the column density contour at 3.2 x lO'^^ 
cm~^ labeled. All of the star formation in M33 is within 
this contour according to the GALEX FUV map (see 
Figure [5]) . Also labeled are a northern arc of gas that 
extends from the northern side of M33's disk and wraps 
back towards the disk and a southern cloud that extends 
in a filament back towards the galaxy. The northern arc 
includes part of the known warp of M33 and extends out 
to projected distances of 11-18 kpc from the contour at 
3.2 X lO'^" cm~-^ 16 — 22 kpc from the galaxy cen- 
ter), while the southern cloud and southern part of the 
warp (to the east of the cloud) extend ~ 10-11 kpc from 
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this contour of the disk. Additional properties of these 
features are summarized in Table 1. Other notable fea- 
tures evident in Figure [T] are the large extent of the HI 
gas (28 X 38 kpc at maximum) and three fingers of gas 
extending away from the extreme edges of the HI distri- 
bution to the east and west. The total HI mass of all of 
the emission shown is Mm = 1.4 x 10^ M©. 

Figure [2] shows M33's HI column density contours over 
the far-UV (FUV; 1350-1750 A) GALEX image (Icr sur- 
face brightness limit of 28.0 mag arcsec"^; Thilker et al. 
2005). The FUV emission is largely limited to the inner 
8 X 15 kpc. The HI thus extends over 2-3 times as far 
as the stellar component traced by the UV light that is 
now commonly used to trace stellar clusters forming in 
the outskirts of galaxies (Werk et al. 2010; Thilker et 
al. 2007). An extension of M33's HI over the optical 
stellar component has also been noted (e.g., Corbelli et 
al. 1989). The FUV emission is confined within the HI 
column density contour at 3.2 x 10^" cm^^ labeled in 
Figure [1] The compact nature of the contours near the 
disk make it difficult to identify which contour the star 
formation begins at, but it can most likely be taken in 
further to the contour at 5.9 x 10^" cm~^. Since the gas 
beyond the 3.2 x 10^° cm~^ contour has a HI mass of 
2.5 X 10^ Mq, over 18% of the HI mass of M33 is beyond 
the star forming disk. 

The complex velocity structure of M33's HI emission 
is best highlighted in the channel maps of Figure [31 
The channel maps have been extended from -412 to - 
46 km (LSR) to encompass a large high- velocity 
cloud known as Wright's cloud (Wright 1974, 1979) and 
emission from the Galactic plane. Each of the ~5 km 
s~^ channels in Figure [3] shows gas with column densi- 
ties as low as 10"'^^ cm~^. The intensity weighted velocity 
of M33's HI emission along a given sightline is also shown 
in Figure[4l All velocities are in the LSR reference frame. 

The first feature to appear in the channel maps of Fig- 
ure [3] is Wright's cloud. Wright's cloud extends from 
the channels at -411 to -339 km s^^ and spatially this 
cloud is known to continue for at least another 1.5° to 
the west of the data shown here (Wright 1979; Lockman 
et al. 2002). The structure of this entire complex as 
mapped by GALEA will be presented in a future paper. 
There is no evidence for a direct link between Wright's 
cloud and M33. They are separated in velocity space by 
~ 20 km s~^, and at this closest point in velocity they 
would be spatially separated by over 50 kpc at the dis- 
tance of M33. In the integrated intensity maps of both 
objects they appear to be separated by only ~ 1.25° (e.g. 
Figurein]), or 16 kpc at the distance of M33, but the ve- 
locity separation at this closest point is over 250 km s~^. 
The mass of the part of the complex shown here is 186 
D2(kpc) M©. 

M33's HI emission begins at a velocity of -324 km 
in Figure [3] and shows neither a diffuse link to 
Wright's cloud nor signs of a filament in the direction 
of M31. The emission subsequently extends both south 
into the main body of M33's disk and north into M33's 
known warp and the northern arc labeled in this figure 
and Figure [TJ A finger of HI emission extends off of the 
northern arc at ^ —241 km s~^ and the arc becomes 
basically indistinguishable from disk emission at about - 
200 km s^^. The mass of the entire northern arc starting 



from where it emanates from the elliptical disk contours 
in Figure[T]is 7.7 x lO'^ Mq. 

Continuing through the channel maps, much of the 
HI emission surrounding the high column density gas in 
M33's disk is difficult to classify into individual features. 
This emission is diffuse in that it surrounds the major- 
ity of the disk and is at a relatively low column density 
compared to the disk emission. In particular, this diffuse 
emission surrounds the eastern side of the galaxy from 
about -200 to -100 km s~^ and eventually merges with 
the southern warp. The southern cloud labeled in this 
figure and Figure [1] is notable as a somewhat more iso- 
lated feature and is most evident in the channels from 
-159 to -138 km s~^. This cloud links back to the disk of 
M33, and the filament from the cloud to the disk can no 
longer be easily isolated by the channel at -118 km s~^. 
This cloud has a HI mass of 1.4 x 10^ M© at 730 kpc 
and w as mapped at high resolution by iWestmeier et all 
()2005[ ). There is another filament in the central regions 
north of the filament to the southern cloud that extends 
into more negative velocities and is obscured by the disk 
emission in the integrated maps. Galactic emission be- 
gins to overlap with M33's emission at velocities less neg- 
ative than -107 km s^^. The last emission which is clearly 
part of M33 is at the tip of the southern warp at -46 km 
s-i. 

Two other figures are included to represent additional 
properties of M33's HI emission. The first is an intensity 
weighted velocity dispersion map of M33 (Figure [5|). The 
gas beyond the star forming disk generally has velocity 
dispersions < 15 km s^^, while the gas in the star form- 
ing disk has a mean velocity dispersion of 18.5 km s"'^. 
The central regions clearly have the highest dispersions 
with values up to 30 km s~^. The decreasing values to- 
wards th e outer regions are foun d for other local gas-rich 
galaxies ([Tamburro et al.ll2009l ). Tamburro et al. assess 
the effect of beam smearing on the velocity dispersion 
of a galaxy (NGC 5055) with a slightly higher inclina- 
tion than M33 (59° vs. M33's 56°) and a 265 pc beam, 
and find the contribution is at most 10% in the outer 
regions and 20% in the inner regions. Finally, Figure [6] 
shows a three dimensional map of M33's HI structure. 
The disk is saturated to bring out M33's surrounding 
features. The velocity range extends from -360 to -52 
km s~^, and therefore includes some Galactic emission 
to the right. The northern arc is clear, as well as the 
southern warp and diffuse gas that connects back to the 
disk. 

3.2. Warped Disk Model Comparison 

We compare our data with the warped disk model of 
CS97 to examine how the features identified here are re- 
lated to this model. The CS97 model is a series of nested 
rings, each with a HI surface density (S (r)), circular ve- 
locity {V {r)), inchnation angle ((/)(r)), and major axis 
position angle (O (r)). These parameters vary with ring 
radius using sigmoidal functions and a total of 14 inde- 
pendent parameters are used to specify these functions as 
outlined by CS97. To directly compare this model with 
our HI data set we also required a line width parameter 
(cr) and a point-spread function (PSE). We use the stan- 
dard Arecibo PSF of 3.4' EWHM. To determine a we ran 
a least-squares fit of the model to the data allowing only 
a to vary, using the MPFIT.PRO program within IDL 
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(Markwardt 2009). We find a best fit a = 14 km s^^, 
very close to the value adopted in CS97, a = 12 km s^^. 

Comparisons between the warped disk model and the 
GALFA-HI data are shown in Figures [7] & [S] As noted 
by CS97, the warped disk (inclined ~ 30° to the central 
disk) reproduces the twisted overall velocity structure of 
M33, but an additional component is needed to accu- 
rately reproduce M33's HI distribution in all areas. The 
model over-predicts the HI flux in the northeast, and to 
a lesser extent in the south. The most significant places 
where the model cannot reproduce the observed HI struc- 
ture are in the majority of the features discussed above, 
including: the northern arc, the southern cloud, the fil- 
aments in the central regions (see Figure [7]), and the 
diffuse emission and fingers to the east and west. There 
is also an offset in the velocity of the model relative to 
the data on the east and west sides of the galaxy. For 
instance, in Figure [8] if one of the six central model con- 
tours is followed from the east to the west, the underlying 
observed HI emission is at a different velocity (color) on 
each side. This may represent the mixing in of the gas 
beyond this model radius. The gas detected beyond the 
model radius tends to have more negative velocities than 
the gas detected within this radius to the west, and more 
positive velocities to the east. 

3.3. Cataloguing Discrete Features 

Isolating halo gas from disk gas is difficult and can 
introduce false separations of related features. Besides 
the visual inspection of features described above, we also 
approached the problem of defining discrete features by 
running an automated object finder called Duchamp-'^^ 
on the data. This allowed us to search for faint discrete 
objects that may have been missed by eye and provided a 
quantitative method of distinguishing objects. Duchamp 
allows one to specify a detection threshold (initial seed 
of a detection must be above this level) and merges fea- 
tures that are linked in position and/or velocity at this 
threshold, or a different threshold if desired. Duchamp 
therefore combines features that are linked by diffuse fil- 
aments and catalogs truly discrete features as individual 
objects. 

When running Duchamp we first applied a banning 
smooth with a width of 5 channels to increase the chances 
of detecting low column density clouds with velocity 
widths similar to Galactic HVCs (~25 km s~^).^^ We 
then ran several iterations of Duchamp and found using 
a detection threshold and merge parameter of 0.08 K (3cr 
for this smoothed cube) was the most successful in terms 
of finding truly discrete clouds and not pushing too far 
into the noise. The results of this run are shown in Fig- 
ure [9] and tabulated in Table 2. Because of the increase 
in the noise level as Galactic emission is approached, we 
cut the search off at —100 km s~^. We therefore do 
not calculate the total flux of M33 (object #11 in Fig- 
ure (5]) with this method and exclude two clouds detected 
at the edge of the velocity range that were clearly part 
of a Galactic filament (see Figure [3]) . We extended the 
search to —450 km s~^ in the other velocity direction to 
include the clouds associated with Wright's cloud. Each 

http:/ /www. atnf.csiro.au/people/Matthew.Whiting/Duchamp/ 
We also ran Duchamp on the original velocity resolution cube, 
but this did not catalog additional discrete objects. 



object catalogued by Duchamp was examined in the data 
cube to assess if it is a genuine detection. This led to the 
removal of objects # 8 and 10 from Table 2 and Figure [9l 

Independent of M33 and Galactic emission, 12 discrete 
real objects were catalogued by Duchamp between -450 
to -100 km s~^ in the region shown in Figure [9] (Table 
2). 7 objects have central LSR velocities between -400 
to -350 km s~^ and are almost certainly associated with 
Wright's cloud (bottom half of Table 2). Of these 7, 
objects #1,3, and 7 are in a noisier part of the cube 
and as can be seen from Figure[3]could be considered part 
of object #2, or the main part of Wright's cloud. Object 
#5 appears spatially removed from Wright's cloud here, 
but is most likely also associated, as an extension in this 
direction was also noted by Lockman et al. (2002) and 
Braun & Thilker (2004). For the remaining 5 objects at 
the top of Table 2, #9 is clearly related to the northern 
arc of M33 and #12 to the southern cloud. The other 3 
objects (#13, 14 & 15) are at velocities between -130 to 
-113 km s^^ and are further from M33 (6-19 kpc from 
M33's HI edge in this figure at the distance of M33). If 
located at the distance of M33 (730 kpc), the objects in 
the first part of Table 2 have an angular extent of 1-3 kpc^ 
velocity widths of 17-30 km s^^, and masses between 10 
and 2 x 10^ Mq. 

Though the data presented here only continue to de- 
clinations of -1-32.5° at full sensitivity, we find no fea- 
tures north of what is shown in Figure [9] beyond emission 
which is part of Galactic filaments. The diffuse filament 
towards M31 claimed by Braun & Thilker (2004) is not 
continuous in velocity with M33 and has LSR velocities 
generally less negative than -100 km s~^. This detec- 
tion may be part of the diffuse filaments extending into 
our Galaxy partially shown in the channel maps of Fig- 
ure [3l The heavy kinematic and spatial smoothing used 
in Braun & Thilker (2004) may have blended parts of 
the Galactic emission into their maps. Future deep map- 
ping of the region between M31 and M33 will resolve the 
origin of any clouds found in this region. 

4. ORBITAL HISTORY OF M33 

To explore the origins of M33's gaseous features we in- 
vestigated the feasibility of a tidal interaction between 
M33 and M31 by constraining the orbital history of the 
system. Given the large present separation between M33 
and M31, a significant tidal interaction that could lead 
to the gaseous features seen here is only possible if past 
orbits brought the two galaxies closer to each other. This 
modeling also puts constraints on the presence of satel- 
lites around M33, as they would not remain bound to 
M33 if it had a previous close encounter with the more 
massive M31. 

The orbital history of M33 can be constrained by in- 
tegrating the motion of M33 backwards in time through 
M31's evolving gravitational potential. Since the proper 
motion of M33 has been recently measured (190 ± 59 km 
s~^ relative to Earth; Brunthaler et al. 2005) and the ra- 
dial velocities of both galaxies are well known (—39 km 
s~i for M33 and -116 km s"^ for M31 relative to the 
Galactic center) , the only unknown velocity components 
are the tangential velocities of M31 in the observed ref- 
erence frame. We calculated a large set of possible orbits 
for the two tangential velocity components of M31, rang- 
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ing each from —200 to +200 km s^"'^. Higher velocities 
are unlikely in a relatively poor group of galaxies such as 
the Local Group. Our grid selection method is similar 
to that in Loeb et al. (2005). We treat these possible 
orbits as a statistical ensemble and use it to estimate the 
probability of a recent close encounter. Following van der 
Marel & Guhathakurta (2008), we use a distance of 770 
kpc for M31 and 794 kpc for M33. The variation in the 
distance for M33, compared to the 730 kpc from Brun- 
thaler et al. (2005), does not affect our results within 
the uncertainties. 

For the purpose of orbit calculations, we treat M33 as 
a point mass moving in the combined gravitational po- 
tential of M31 and the Galaxy. The mass distribution 
of M31 is modeled as the sum of a baryon disk and a 
dark matter halo, with the total virial mass Mvir(O) = 
2 X 10^^ Mq at redshift 2 = 0. We adopt the disk mass 
Md{0) = 4 X 10^° Mq. The disk is described by an ana- 
lytical Miyamoto-Nagai profile with a scale length Vd — 5 
kpc. Although our adopted baryon mass may be an un- 
derestimate of the M31 disk and bulge, we show below 
that it plays a much smaller role than the dark halo in 
determining the possible orbits of M33. A more detailed 
model of the M31 stellar distribution by Fardal et al. 

(2007) had a similar scale length but larger total baryon 
mass. 

For the halo we adopt a NFW profile with a concen- 
tration parameter c=12atz = 0, in line with the 
expectations of ACDM simulations. This results in a 
scale radius Vs ~ 25 kpc. The halo mass distribution 
is the main uncertainty of our calculations. Unfortu- 
nately, the outer mass distribution of M31 is not well 
determined observationally. Fardal et al. (2007) and 
Seigar et al. (2008) calculate Mvir « lO^^ Mq, while 
Loeb et al. (2005) and van der Marel & Guhathakurta 

(2008) adopt Afvir w 3 x 10^^ Mq. Such widely different 
estimates can be partly due to the different analysis and 
partly due to the different assumed profiles. In particu- 
lar, Fardal et al. (2007) find the best fit to the rotation 
curve requires the concentration parameter c ~ 30, much 
higher than that expected for AGDM halos of such mass. 
A very concentrated lower-mass halo and a less concen- 
trated higher-mass halo can contain similar amounts of 
mass within a given radius. For example, in our model 
the total masses of stars and dark matter enclosed within 
60, 120, and 200 kpc are (0.66, 1.15, 1.60) x lO^^ jyj^^ re- 
spectively, while in the Fardal et al. (2007) model these 
masses are (0.51,0.72,0.88) x 10^^ ^s we find no 
plausible orbits of M33 coming closer to M31 than 60 
kpc, the size of the baryon distribution in unimportant 
for this calculation. 

Since the mass of a galaxy is expected to grow with 
time due to hierarchical merging and accretion, we ac- 
count for this evolution by reducing the mass as we inte- 
grate the orbits in the past for 10 Gyr from the present 
time, to. We assume that the disk mass scales linearly 
with cosmic time, Md{t) — Md(Q)t/tn. For the halo 
mass , we adopt an exponential growth (jWechsler et al.l 
I2002D with the formation redshift z ss 1, which gives 
Mvir{z) = Mvir(O) e~^. Note that for most orbits we find 
the last close passage between M33 and M31 within the 
last 3 Gyr {z < 0.3). In this relatively short interval the 
gravitational potential is not expected to have changed 



significantly, and therefore, our assumptions about its 
evolution with redshift are not affecting our conclusions. 

In addition to M31, we include the gravitational po- 
tential of the Galaxy, treating it as a point-mass (10^^ 
Mq) with an analogous redshift evolution. It turns out 
to have a negligible effect due to the large current dis- 
tance to M33 and an even larger separation in the past. 
We also include dynamica l friction on the orbit of M33 
(jBinnev fc Tremaiiiil2008D . but the effect is again small 
because the galaxy never gets close enough to M31. 

The main goal of the orbit calculations was to deter- 
mine the evolution of the tidal radius of M33, r^. Any 
material outside rt is no longer gravitationally bound 
to the small galaxy and would be pulled away by tidal 
forces. However, it should be noted that the isopotential 
surfaces are not spherical and a single value of the tidal 
radius gives only a rough idea of the tidal truncation 
of M33. To estimate rt we assume flat rotation curves 
for both M33 and M31 around rt, with the circular ve- 
locities Vm33 ~ 100 km s~^ and Vmsi w 250 km s~^, 
respectively. This gives vt w 2/3^/^ i?p (Vmss/Vmsi) ~ 
0.15 i?p, where Rp is the perigalactic distance between 
M33 and M31. 

The results of the orbit calculations are shown in Fig- 
ures [TUHI31 Figure [TU] shows the probability (cumula- 
tive fraction of orbits) of a given tidal radius for the 
orbits with the maximum tangential velocity of M31 re- 
stricted to be below 50, 100, and 200 km s^^. Although 
about 30% of the simulated orbits do not have a pre- 
vious closer approach, these orbits seem unlikely as it 
would imply that M33 did not form in the vicinity of 
the other galaxies of the Local Group. Figure [10] shows 
that at all possible M31 velocities leading to a previous 
closer approach between M33 and M31, there is a 60% 
probability of the galaxies once being within 100 kpc of 
each other, resulting in a tidal radius of M33 below 15 
kpc. For the most likely scenario of V(M31)ta„ < 150 
km s^^ (see also Ivan der Marel fc Guhathakurtal ()2008f ) 
who find V(M31)ta„ < 56.3 km s'^), there is more than 
a 75% chance of the M33 tidal radius being smaller than 
15 kpc. Even if we include the orbits lacking a previous 
closer approach, at these most likely M31 velocities the 
probability for rt < 15 kpc is still very high, about 60%. 
Figure [TT] shows a 3D representation of one of such orbits 
with rt « 10 kpc. 

Figure [T2l shows the probability distribution of the dis- 
tance and lookback time of the close encounters for the 
entire range of velocities. The most likely encounters 
are the closest and most recent (darker shading on the 
plot). Approximately a third of all possible orbits lead 
to Rp < 80 kpc and rt < 12 kpc, which occurred between 
1 and 2.5 Gyr ago, and half result in a close encounter 
resulting in rt < 15 kpc within the last 3 Gyr. Even 
though this is a statistical argument, it provides support 
for the idea of a relatively recent tidal interaction be- 
tween M33 and M31 even if the M31 velocity is much 
higher than expected. Figure [T2] shows the probability 
distribution of the M33 velocity relative to M31 at the 
last close encounter. The majority of the orbits indicate 
a velocity around 400 km s~^. 

5. DISCUSSION 

M33's gaseous features are connected to the galaxy 
in position and velocity and cause M33 to extend over 
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28 X 38 kpc at maximum extent. Given the likelihood 
of a previous closer approach with M31 resulting in a 
tidal radius < 15 kpc (see Section 4), these features were 
most likely pulled from the outer disk of M33 during this 
passage. The gaseous disks of spiral and dwarf galaxies 
are commonly found to extend significantly further than 
the st ellar components ()Meurer et al.lll996t iBegum et al.l 
|2005[) and will be the first thing stripped from a galaxy 
during an interaction. The early discovery of M33's warp 
was an initial hint of M33's disruption (Rogstad, Wright 
& Lockhart 1976; CorbeUi & Schneider 1997), but at that 
time it was concluded that the tidal force of M31 rela- 
tive to M33 was not strong enough to cause the warp. 
The recent proper motion measurements of M33 have 
permitted the orbit of this galaxy relative to M31 to be 
assessed. Our analysis, as well as the analysis of Loeb et 
al. (2005) and van der Marel & Guhathakurta (2008), 
indicate a previous interaction between M33 and M31 is 
highly likely. The tidal radius of < 15 kpc and approach 
of < 100 kpc within the past 3 Gyr occurs in more than 
75% of the most likely orbits (Vm31 < 150 km s~^) in 
our analysis. Given M33's gaseous features extend out 
beyond this radius and the largest features are relatively 
symmetric, a tidal origin seems likely. In addition, since 
a satellite would not remain bound to M33 unless located 
at < 15 kpc, the gaseous features are unlikely to be the 
remnants of an accreting satellite. Deep stellar surveys 
have also not revealed any clearly distinct stellar features 
in the vicinity of M33 (e.g., Ibata et al. 2007). 

As M33 passed through M31's extended halo at most 
likely distances of 65-120 kpc (Figure [T^]), ram pres- 
sure forces from M31's hot gaseous halo are also likely 
to have acted on M33's gas. Though there is no di- 
rect evidence for a hot halo around M31, it is inferred 
to exist thro ugh the la ck of gas in the dwarf satel- 
lites of M31 (iGrcevich fc Putmaiil [20091 ) . head-tail halo 
clouds (jWestmeier et al.i .2005.) , and the predictions of 
simulations of galaxy formation (jSommer-Larsenl 120061 : 
iKaufmann et al.''2006'). Using the simple Gunn & Gott 
criterion (|Gunn fc Gott 1972) as used by VoUmer et al. 
(2008), 
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we can estimate if M31's hot halo is likely to have played 
a role in stripping the outer parts of M33. M33 would be 
moving at approximately Vgai — 300 — 415 km s~^ rel- 
ative to M31 at perigalacticon fFigure [T3|) . Using a ro- 
tatio n al velocity of Vrot = 100 km s~^ (jCiardullo et al.l 
12004 iCorbelli fc Saluccil l2000l) . a stripping radius (R) 
of 9 kpc (lower bound corresponding to the edge of the 
FUV disk), and a T,gas of 3 x 10^° cm'^ (see Table 1) 
suggests M31 halo densities greater than 10~'^'^ cm~'^ are 
needed for ram pressure stripping to be effective on the 
outer regions of M33. At distances of 65-120 kpc a hot 
halo around galaxies similar t o M31 is highly unlikely to 
have these types of densities (Bregmanl l2007t iPeek et aD 
I2007t [Gaensler et al. 2008). Densities on the order of 
10~^'^ — 10~'^'^ cm~'^ are expecte d from the galaxy for- 
mation simulations at these radii (|Sommer-Larsenl|2006l : 
IKaufmann et al.ll2006( ) and the observational constraints 
referenced above are most consistent with the lowest 
end of this range. Therefore, the role of a diffuse halo 



medium around M31 in stripping out M33's gaseous fea- 
tures is likely to be minor, consistent with their symmet- 
ric structu re not be ing indicativ e of ra m pressure strip- 
ping ( VoU mer et al.! 2008; Chung [2007[). 

Though ram pressure effects are unlikely to lead to the 
origin of M33's gaseous features, they will apply a drag 
on this relatively low column density gas (compared to 
the disk) at large radii from the center of M33's gravi- 
tational potential. The fingers extending from the sides 
of the galaxy may be evidence of this effect. The combi- 
nation of this drag and the tidal forces will cause M33's 
outer gas to lose angular momentum and fall towards the 
central regions of M33 as it moves away from the gravita- 
tional influence of M31 (iMoore fc pavislll994t iGardineii 
fT99llBarnesl[200llCox et al.ll2008n . None of the HI gas 
surrounding M33's stellar disk is above M33's escape ve- 
locity and its infall may be subtly evident in its kine- 
matics. The gas on the western side of M33 has more 
negative velocities than the gas in the underlying disk, 
consistent with an infall of gas from the far side of the 
galaxy; while the gas on the eastern side has more posi- 
tive velocities than the gas in the disk consistent with a 
net inflow from the near side of the galaxy. Gas falling 
back onto M33 may represent an important method of 
fueling and triggering the star formation in M33's disk 
(see Section 5.2.2). 

5.1. M33 Halo Clouds 

The limited number of discrete objects catalogued by 
Duchamp in the vicinity of M33 (Section 3.1) high- 
lights that the majority of the surrounding gas is di- 
rectly connected to M33's disk in position- velocity space. 
This include s the majority of the clouds catalogued by 
IGrossi et al.l ((2008). The clouds that are catalogued as 
discrete objects (see Figure [5]) are either very close to 
being connected to M33 in position-velocity space (#9, 
12), associated with the Wright cloud complex, or are the 
three clouds at low velocities (Vgsr = 4 — 20 km s~-^; 
#13, 14, 15). These lower velocity clouds may be part 
of a Galactic filament or they may be related to M33, as 
they do extend towards the position and velocity of a fil- 
ament on that side of the galaxy. If related they extend 
~19 kpc from the HI edge of M33 (~38 kpc from the 
center) and may represent additional material disrupted 
during the interaction with M31, or possibly some type 
of cold flow of gas towards the disk (Keres et al. 2005). 
Deeper observations will help to confirm their relation- 
ship, if any, to M33. 

The properties of the clouds that are associated (#9, 
12) or potentially associated (#13,14,15) with M33 (top 
part of Table 2) can be compared to the halo clouds of 
our Galaxy and M31 (referred to as high- velocity clouds 
(HVCs) for our Galaxy). The clouds have similar prop- 
erties to HVCs and M31 halo clouds in terms of their col- 
umn densities and velocity widths (iPutman et al.ll2003 
Ide Heij et all[200l lm"stmeier et al.ll2005n . The similar 
velocity width of the clouds (20-30 km s""'^) is indica- 
tive of a warm cloud component (~10,000 K) and sug- 
gests a simil ar nature for the M33, M 3 1 and Milky Way 
halo clouds (iStanim irovic e t al.l [20081 iBen Bekhti et all 
[200l iWestmeier e~ al. 2005). For Galactic HVCs this 
warm component is often interpreted as representative 
of a pressure equilibrium with a surrounding hot halo 
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medium (jWolfire et al.lllQQSf ). M33 is not generally ex- 
pected to have a remnant primordial hot halo as it is 
not considered massive enough for a large quantity of 
gas to be shock- heated as it flows into the dark mat - 
ter potential well (Keres et al1l2005l : iBrooks et al.ll2009n . 
In addition, if there was any surrounding halo medium 
it would most likely have been stripped during the en- 
counter with M31. Therefore the warm cloud compo- 
nent may represent warm material stripped from M33, 
or heating from either the interacti on or esc aping ra- 
diation from M33's disk (Bland-H awthorn fc Maloney 
11999') . Sensitive Fabry-Perot Ha observations of the ex- 
tended gaseous M33 system would be a very interesting 
complement to the data presented here. 

At the distance of M33, the discrete clouds also have 
similar masses to the M31 clouds and the l arge Galactic 
HVCs with distance c onstra i nts (~ 10 kpc; iThom et al.l 
(I2OOI [2008h : [Wakker eFall (pOOSh ): but the physical 
sizes tend to be much smaller. This may be consistent 
with the largest clouds catalogued as halo clouds in M31 
and the Galaxy either showing a connection to the disk, 
or in the case of the Galaxy, the fact that they would 
show a connection if observed by someone in an external 
galaxy (jThilker et alllmillKaiijerla et al.|[200l . For in- 
stance, for M33 the southern cloud is connected to the 
gaseous disk (see Figured]), and is more typical in size 
(though still on the small side) to Galactic HVCs. Its 
mass would be similar to some of the larger Galactic 
HVC complexes. In any case, the clouds surrounding 
M33 show numerous similar properties to the Galactic 
and M31 halo clouds and this could be suggestive of a 
similar origin or destruction mechanism. 

The clouds catalogued as part of the Wright's cloud 
complex are unlikely to be associated with M33 given 
their projected distance from the galaxy, velocity off- 
set, and the large mass and siz e of the compl ex at the 
distance of M33 (> 10« M©; IWrightl dHzl)). This 
cloud may be associated with the tail of the Magellanic 
Stream, as the Stream is found to extend quite far in 
the direction of M33 and has sever al features at simi- 
lar extreme negative LS R velocities (jStanimirovic et aTl 
120081: IWestmeier fc Koribalski ,200a) . At 60 kpc this 
cloud would have a mass on the order of 10^ M©. The 
linewidths of the catalogued clouds that are part of the 
Wright's cloud complex are again typical of Galactic 
HVCs, as is the angular size of the primary part of the 
complex. The GALFA data of this complex will be ana- 
lyzed further in a future paper with a similar method to 
that outhned in Peek et al. (2007). 

5.2. Future of M33 
5.2.1. Fueling M31 

The most likely fate of M33 is to fuel M31 given the 
current disruption of this galaxy and their likely mo- 
tion relative to each other. M33 will therefore provide 
1.4 X 10^ M©-^^ in star formation fuel to M31, which cur- 
rently has only 5 — 6 x 10^ M© in fuel. Depending on the 
orbits of the galaxies, M33 may eventually resemble the 
Magellanic System, with leading and trailing streams of 
gas being stripped through some combination of tidal and 

The molecular component of M33 is estimated to be only an 
additional 4.5 X lO'' Mq (lEngargiola et al.„200a) . 



ram pressure forces (iPutman et al. l l2003t iConnors et al.l 
l2006t iMastropietro et al.ll2005f ). 7t is possible that some 
of the clouds in M31's halo were left behind from the first 
passage of M33. In particular there are several clouds 
at projected distances up to 50 kpc from M31's disk 
that c annot be associated with any known stel lar fea- 
tures (IWestmeier et al.ll2Q07l: pThilker et al]|200l . 

Mergers with large satellites may be an important 
source of star formation fuel for spiral galaxies like M31 
and the Milky Way. The Magellanic System and the 
Milky Way are even closer in total gas mass than M33 
and M31, with the Magellanic System host ing ~ 10^ M©, 
close to 50% of the Milk y Way's gas mass (jPutman et al.l 
l2003t iLevine et al.l|200 6l . Chemical evolution models for 
the Gal axy typically require an average inflow of ~ 1 
M©/yr (jChiappini et all l200ll : iFenner fc Gibsolll [200l . 
and obtaining all of this from lower mass satellites, or 
gas condens i ng fro m the hot halo a ppears to be difficult 
(iPeek et al.l 120081: iPutmanI [200l iGrcevich fc PutmanI 
|2009() . The equivalent of one system of similar gas mass 
to the Magellanic System or M33 per Gyr over the past 
5-7 Gyr would have provided a galaxy like the Milky Way 
or M31 with enough fuel to maintain its star formation 
and not become severely depleted. 

A scenario of larger, gas-rich, satellite mergers pro- 
viding the majority of the star formation fuel for a 
spiral galaxy may also be consistent with the sur- 
vival of a galaxy disk according to recent simulations 
(Robertso n et all 120061 : IBrook et all 120041 : Poster et all 
,2009) . Though the M31:M33 total mass ratio is not 
quite at the 1:10 mass ratio often considered particu- 
larly harmful to the disk, at approximately 1:20-40 it 
is significantly larger than the majority of the dwarf 
satellites generally considered as fuel for Local Group 
galaxies. The simulations generally find that when the 
merging system is gas dominated, the thin disk is ei- 
ther able to re-form as the gas settles back onto the disk 
and/or the kinetic impact energy is absorbed by the gas. 
The evolution of satellites towards radial orbits as they 
merge may also sub stantially suppress the disk heating 
(jHopkins et al.l l2008) , although the addition of gas seems 
to be a more widel y acc epted result in aiding disk sur- 
vival (jPurcell et ahfeoOQl ). The presence of M33 and the 
Magellanic System in the context of M31 and the Milky 
Way's evolution indicates these types of gas- rich mergers 
continue to occur at z = 0. 

5.2.2. Fueling M33 

M33 is currently form ing stars at ^ 0.7 M©/yr 
(jBlitz fc Rosolowskvl[2006l) . The galaxy's HI mass indi- 
cates it only has 2 Gyr left of star formation activity at 
this rate unless a source of fuel is obtained from the Local 
Group. Thus M33 appears to be in a similar situation to 
the Milky Way and other local spiral galaxies (e.g.. Peek 
2009); it is either in a unique place in time where it will 
soon run out of star formation fuel, or it has a continuous 
mode of obtaining star formation fuel from a surrounding 
medium. Recent chemical evolution models of M33 also 
suggest continuous gaseous inflow is needed over the life- 
time of the galaxy to expla i n the range of stellar meta l- 
licities (jMagrini et al.ll2007t [Barker fc Saraiedinill2008[ ). 

A source of continuous fuel for M33 is difficult to find. 
M33 may be partially fueling itself with the gaseous fea- 
tures that extend beyond the current star forming disk. 
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but this is only 2.5 x 10^ Mq and would rapidly be ex- 
hausted. Satellites cannot provide M33 with fuel given 
the close passage with M31 would have stripped them 
from MBS. As mentioned previously, the lack of satellite 
accretion for M33 is also consistent with its smoot h stel- 
lar d istribution in deep surveys to date (Ferguso n et al.l 
I2007D . The replenishment of M33's HI through cool- 
ing of gas near the disk is a possibility given the con- 
nection of M33's gaseous features to the star forming 
disk. The Milky Way, M31, and other spiral galax- 
ies also tend to have the majority of their halo gas 
close to the disk, and the higher densities there make 
cooling of a diffuse halo medium to HI a more likely 
possibility (iHeitsch &: PutmanI 120091 : iPeek et al.l 120081 : 
iFraternali &: Binnevll2008[ ). If this is the case, the source 
of the cooling gas remains a mystery, as the passage with 
M31 is likely to have stripped any M33 halo medium and 
there is no strong evidence for a substantial diffuse Local 
Group medium. 

Despite the limited mass present in M33's gaseous 
features for long term feeding, on the short term the 
inflow of this gas back towards M33 may be affect- 
ing its star formation. The pressure of the infalling 
gas may help to explain the efficiency of M33's cur- 
rent star formation. Several groups have found that 
M33 is more efficient at converti ng molecular gas to 
stars than typical nearby spirals (Gar dan et al.l [2"007l : 
Blitz & Rosolowskv 2006; Kcnnicutt 1998). For instance, 
Blitz fc Rosolowskyi (|20061 ) are able to predict the molec- 
ular gas fraction with their pressure regulated star forma- 
tion prescription, but under-predict the star formation 
rate by a factor of 7. An additional pressure on the gas 
generated by material falling back onto the disk and the 
higher velocity dispersion in the disk than assumed for 
all galaxies by Blitz & Rosolowsky (2006) may be a com- 
ponent in the discrepancy. Blitz & Rosolowsky assumed 
a constant velocity dispersion of 8 km s~^ across the 
disks of all of their galaxies, and their mid-plane pres- 
sure is directly proportional to this value. Still M33's 
higher mean value of 18.5 km (see Figure [5|) should 
affect the molecular gas fraction as well as the SFR un- 
der their criterion, so this is not a simple solution to the 
discrepancy in SFR. We note that M33's mean velocity 
dispersion is also elevat ed compared to the sani ple of gas- 
rich galaxies studied bv lTamburro et ahl (|2009D . with the 
exception of NGC 5194, which they argue is elevated due 
to the enhanced star formation triggered by the tidal in- 
teraction with NGC 5195. Since M33 is now beyond the 
dominant gravitational influence of M31, perhaps it will 
soon cycle out of an elevated star formation phase. 

M33 does not appear to show star formation in its sur- 
rounding gaseous structures, or an extended UV (XUV) 
disk as seen for many galaxies with large gaseous fea- 
tures beyond the main optical compone nt (see Figure [2j 
IWerk et al.l (|20Toh : iThilker et al.l (l2007h ). We compared 
our observed column densities at the largest radii M33 is 
forming stars to the expected critical surface density de- 
rived by Kcnnicutt (1998) based on the Toomre (1964) 
criterion and confirmed previous findings that the ob- 
served HI column (3-5 M0 pc~^) is three times lower 
than this expected critical density. Therefore M33 is 
forming stars in relatively low column density gas, but 
for some reason this star formation does not continue 
into the extended gaseous features like it does in other 



systems. The presence (or lack of) star formation in ex- 
tended gaseous features will be an ongoing area of study, 
with M33 being an important local probe of this process. 

6. CONCLUSIONS 

We have presented new HI data of M33 and its sur- 
roundings from the GALFA-HI Survey and have high- 
lighted several distinct gaseous features that extend out 
to 18 kpc from the edge of the stellar component traced 
in FUV emission (Figure [T] & [21 Section 3.1). The max- 
imum extent of M33's HI is 28 x 38 kpc, and includes 
features such as a northern arc of gas that extends from 
M33's warp back towards the disk, a southern cloud with 
a filament to the galaxy, and diffuse emission and several 
fingers to the east and west of the galaxy (Figures 1-6). 
The HI gas beyond the star-forming disk contains 18% of 
M33's total HI mass of 1.4 x 10^ Mq and the vast major- 
ity of the gas surrounding M33 is directly connected to 
the disk in position-velocity space. The GALFA-HI data 
is compared to the warped disk model of CS97 to high- 
light the relationship of the extended gas to this model 
(Figures El & [3 Section 3.2). The following additional 
results are described in the paper. 

• The orbit of M33 relative to M31 is examined to 
investigate the feasibility of a tidal interaction be- 
tween the two galaxies resulting in M33's gaseous 
features (Section[4]). Consistent with previous anal- 
yses of the likely orbits, we find a previous close 
encounter (<100 kpc) between the two galaxies is 
highly likely (> 75% of the orbits for the most likely 
V(M31)tan < 150 km s~^), resulting in a tidal ra- 
dius for M33 of < 15 kpc. This suggests a satel- 
lite of M33 would not remain bound to it and the 
gaseous features were most likely tidally pulled di- 
rectly from M33 in the past 3 Gyr. 

• The combination of tidal and ram pressure forces 
acting on M33 as it passes through an extended 
halo of M31, is likely to cause the tidally extracted 
gas to lose angular momentum and fall back to- 
wards M33's disk as it moves away from M31. This 
is consistent with the velocity sense of the observed 
gas and the elevated star formation in M33's disk. 

• An automated object finder is run on the M33 cube 
to define truly discrete objects from those directly 
connected to the galaxy (Sections 3.3 & 5.1). Be- 
sides the clouds associated with the Wright's cloud 
complex (which is unlikely to be associated with 
M33 given its offset and large size at M33's dis- 
tance), there are only 5 small clouds catalogued, 
and all other gas is connected to M33. Of the 5 
clouds, 2 are very close to being connected to M33, 
and the other 3 are either in a filament going out 
38 kpc from the center of M33 or part of a Galactic 
filament given their lower velocities. These clouds 
show similar properties to Galactic HVCs and M31 
halo clouds, including their velocity widths (^ 25 
km s~^) which are indicative of a warm cloud com- 
ponent. 

• M33 will eventually become star formation fuel for 
M31, and if it does not substantially deplete its gas 
content before then, it will provide ~25% of M31's 
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current HI mass. This together with the Magel- 
lanic System, which hosts close to 50% of the Milky 
Way's current HI mass, indicates the contribution 
of gas from large mergers can play an important 
role in the ongoing fueling of galaxies. The fact 
that satellite mergers that are gas-rich are not ex- 
pected to be overly disruptive to the disk is also 
support for these types of mergers providing sub- 
stantial fuel for a spiral galaxy. 

• M33 will either run out of star formation fuel in the 
next 2 Gyr if it continues to form stars at its current 
rate, or it must accrete gas from an unknown Local 
Group medium. Other notable characteristics of 
M33's gas and star formation are the high mean 
velocity dispersion across the disk 18.5 km s^^; 
Figure[5l) that may be consistent with the relatively 
high SFR and recent interaction, and the fact that 
M33 does not have star formation as traced in the 
FUV in its extended gaseous component, or beyond 
the contours at 3 — 6 x 10^" cm~^ (Figure [21). 
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TABLE 1 
Properties of M33's HI 



Region 


RA=» 


Dec=» 


Vlsr" 


HI Mass ^ 


Nhi (peak) 








km s ^ 




_2 

cm 


Entire Galaxy 


01:33:42 


30:35:30 


-147 


1.4 X 10^ 


1.9 X 10^1 


Gas Beyond Star Forming Disk'^ 








2.5 X 10** 


3.1 X 10^0 


Northern Arc*^ 


01:31:50 


31:18:30 


-232 


7.7 X 10^ 


3.1 X 10^0 


Southern Cloud" 


01:32:30 


29:35:30 


-158 


1.4 X 10'^ 


3.7 X 10^^ 



^ Approximate central value. 

^ At the distance of 730 kpc. The dependency is Mhi ck D^, so if the largest M33 distance of 964 

kpc (|B onanos ct al.ll2006l) is adopted all HI masses should be multiplied by 1.7. 

^ Gas beyond the ~ 3.2 X 10^" cm~^ contour in Fig.[T] 

^ The full arc as it emanates from the regular disk structure in Fig. [T] 

^ Only the discrete part of the cloud evident in Fig. i.e., does not include the filament back to 
the central regions of the disk evident in Fig. [3] 
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TABLE 2 

Discrete HI Sources at Vlsr < —100 km s~^ in Figure[9] 



Obj # 


RA 


DEC 




^DEC 


ARA^' 


ADec'* 




VLSi? 






jTAv 










1 


1 


kpc 


kpc 


K 


km 


km 


km 


K km s-l 


M© 


9 


01:28:42.8 


+31:41:23.5 


7.7 


7.0 


1.63 


1.49 


0.25 


-264 


-123 


30.0 


237 


1.6E+05 


12 


01:33:25.3 


+29:26:50.1 


6.1 


6.0 


1.30 


1.27 


0.22 


-150 


-16 


29.7 


148 


9.7E+04 


13 


01:39:33.1 


+30:41:21.0 


12.0 


8.0 


2.56 


1.70 


0.20 


-129 


4 


28.5 


202 


1.3E+05 


14 


01:44:01.0 


+31:15:52.4 


9.4 


5.0 


2.00 


1.06 


0.22 


-127 


5 


20.9 


92.5 


6.1E+Q4 


15 


01:40:08.1 


+30:49:17.4 


5.2 


5.0 


1.09 


1.06 


0.24 


-113 


20 


17.3 


55.7 


3.7E+04 


Wright 


Cloud 


Objects 






















1 


01:15:18.5 


+30:49:26.8 


17.2 


11.0 


3.65 


2.34 


0.19 


-391 


-246 


37.6 


236 


1.6E+05 


2 


01:20:21.5 


+29:42:17.1 


182.4 


92.0 


38.73 


19.54 


4.33 


-390 


-249 


17.7 


1.35E+05 


8.9E+07 


3" 


01:14:17.6 


+30:40:46.4 


9.5 


6.0 


2.01 


1.27 


0.17 


-383 


-238 


16.1 


29.8 


2.0E+04 


4 


01:19:02.0 


+31:13:01.5 


13.7 


13.0 


2.90 


2.76 


0.39 


-380 


-236 


24.0 


617 


4.1E+05 


5 


01:36:59.4 


+30:03:25.2 


7.8 


10.0 


1.65 


2.12 


0.18 


-375 


-241 


20.1 


71.5 


4.7E+04 


6 


01:18:02.8 


+28:29:37.2 


108.1 


38.0 


22.96 


8.07 


1.85 


-370 


-230 


17.3 


1.42E+04 


9.3E+06 


7 


01:16:17.9 


+30:37:41.0 


34.4 


19.0 


7.31 


4.03 


0.30 


-383 


-239 


24.9 


816 


5.4E+05 



^ At a distance of 730 kpc. 

^ Located at the edge of the cube. 
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Right Ascension (J3000) 

Fig. 2. — The Far-UV GALEX image of M33 with contours representing HI column density overlaid. The contours are the same as 
Figure[T] The FUV emission is largely confined within the 5.9 X 10'^'' cm"'^ contour (the 12th contour counting from the outside). 
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Fig. 3. Channel maps every 5 km from -411 to -46 km s^^ showing the HI emission from Wright's Cloud, M33, and some Galactic 
emission. The color bar at the right shows the intensity in Kelvin of the greyscale and the velocity is labeled in the upper left. Outside of 
the channels that include Galactic emission, the 3a detection level is 0.1 K, or 10^* cm~^ per 5 km s"'^ channel. 
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Fig. 3.— (cont.) 
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Fig. 3.— (cont.) 
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Fig. 4.— 



The intensity weighted velocity (LSR) of the M33 gas with HI column density contours from Figure [T] overlaid. 
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Fig. 5. — The intensity weighted velocity dispersion map of M33 with HI column density contours from Figure [T] overlaid. Some of the 
scanning artifacts from the drift scans are more evident in this map. 




Fig. 7. — Channel maps showing the warped disk model in contours over the HI data incrementing every 15 km s ^. The grey scale for 
the HI data extends from - 40 K on a log scale and the contours for the model are 0.2, 0.4, 0.8, 1.6, and 3.2 K. 
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32°00' - 




-100 



-150 



in 



-200 ^ 



-250 



29°30 



300 



36™ 34™ 32™ 30' 

Right Ascension (J2000) 



Fig. 8. — The intensity weighted velocity of M33's HI (see Figure |4]l with the velocity contours of the warped disk model overlaid. The 
contours start at -270 km s~^at the top and extend to -75 km s~^at the bottom incrementing in 15 km s~^intervals. 
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distance of closest approach (kpc) 




tidal radius (kpc) 



Fig. 10. — Cumulative fraction of orbits that result in the tidal radius for M33 noted on the bottom x-axis and and the distance of closest 
approach between M31 and M33 (perigalacticon) noted on the top x-axis. The three different pairs of lines show subsets of orbits with the 
maximum tangential velocity of M31 restricted to be below 50, 100, and 200 km s~^, respectively. The solid lines are for orbits that had a 
perigalacticon closer than M33's current separation (~ 200 kpc), and the dashed lines represent the total grid of orbits. 
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Flc;. 11. — A three-dimensional representation of a typical possible orbit of M33 shown in the M31 reference frame. The solid line shows 
M33's orbit when the Milky Way potential is included and the dashed line shows the orbit when it is excluded, 'peri' marks the closest 
approach of M33 to M31 (~ 10 kpc) on this orbit 1.6 Gyr ago. The Milky Way moves outside the plotted box after only 0.5 Gyr in the 
past and therefore has little influence on the orbit. 
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60 80 100 120 140 

Rp(kpc) 

Fig. 12. — The distribution of pcrigalactic distances and lookback times of the closest approach between M33 and M31, for the full range 
of M31 tangential velocities up to 200 km s~^. The contours represent the cumulative percentage of calculated orbits within the noted 
range of perigalactica and lookback times. The contours have abrupt left and bottom sides because there are no orbits with Rp < 63 kpc 
and tp < 1 Gyr. 
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Fig. 13. — The distribution of velocities and lookback times at the closest approach between M33 and M31. The contours represent the 
cumulative percentage of calculated orbits within the noted range of velocities and lookback times. 



